Easy implementable algorithm for the geometric measure of entanglement 
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We present an easy implementable algorithm for approximating the geometric measure of entan- 
glement from above. The algorithm can be applied to any multipartite mixed state. It involves only 
the solution of an eigenproblem and finding a singular value decomposition, no further numerical 
techniques are needed. To provide examples, the algorithm was applied to the isotropic states of 3 
qubits and the 3-qubit XX model with external magnetic field. 
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I. INTRODUCTION 

Quantum entanglement as a fascinating nonclassical fea- 
ture has attracted attention since the early days of the 
quantum theory [T] |2] • In the last decades its importance 
for the quantum information theory has been recognized, 
since entanglement plays a crucial role in almost every 
quantum computational task 

A bipartite pure state is said to be entangled, if it cannot 
be written in the product form 



\i^i^) = \r)®\r 



(1) 



States which are not entangled are called separable. In 
general, the number of parties is n > 2, and fully sepa- 
rable pure states become 



1^: 



scp/ 



(2) 



The theory of entanglement has also been extended to the 
case, when the quantum state is not pure [HIS]- Then a 
mixed state pscp is called separable, if it can be written 
as a convex combination of separable pure states: 



(3) 



with nonnegative probabilities pi, X^i-P* ~ 1- Quantifi- 
cation of entanglement is one of the main research areas 
in quantum information theory [S]. For bipartite pure 
states, the entanglement is usually quantified using the 
von Neumann entropy of the reduced state: 



E{\^^^))=-Tr[p^\og,p^] 



(4) 



= Tib . For multipartite sys- 

tems and mixed states many different measures of entan- 
glement were proposed 015]. In general, a measure of 
entanglement is any continuous function E on the space 
of mixed states p which satisfies at least the following 
properties [5]: 

• E is nonnegative and zero if and only if the state 
is separable. 



• E does not increase under local operations and clas- 
sical communication: 

E{K{p))<E{p), 

where A is any LOCC operation. 

For bipartite mixed states, an important measure of en- 
tanglement is the entanglement of formation Ef. For 
pure states it is defined as the von Neumann entropy of 
the reduced state as given in Q . The extension to mixed 
states is done via the convex roof construction [HIH]: 



Ef (p) = mm^p.Eili^i)) , 



(5) 



where the minimum is taken over all pure state decom- 
positions of p. 

In this paper we will consider the geometric measure of 
entanglement. For pure states it is defined as follows [H]: 



i?G(|^))-l-max|(^|0)|^ 
|0)es 



(6) 



where the maximization is done over the set of separable 
states S. For mixed states p the geometric measure of 
entanglement was originally defined via the convex roof 
construction, in the same way as it was done for the 
entanglement of formation [S]: 



Eg (p) = min "^piEc 



(7) 



with minimization over all pure state decompositions of 
p. Similar measures of entanglement were also considered 
earlier in (TDl HI] . 

If p is a two-qubit state, general expressions for Ef and 
Eg are known [iinilB]: 



Efip) = h[\ + yi-C{pf 



EG{p)^\[l-\/l-C{py 



(8) 
(9) 
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The concurrence C (p) is given by 

C (p) = max {0, Ai - A2 - A3 - A4} . 



(10) 
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where Xi are the square roots of the eigenvalues of 
p • p in decreasing order, and p is defined as p = 
{ay (g) ay) p* {ay (g) ay). 

For most quantum states no exact expression for any 
measure of entanglement is known, and thus numeri- 
cal algorithms must be used. One of the first algo- 
rithms computing entanglement has been presented in 
[14| . There the entanglement of formation was approx- 
imated using a random walk algorithm on the space of 
the decompositions of the given mixed state. A much 
faster algorithm for the entanglement of formation was 
presented in [TS]. This algorithm made use of the conju- 
gate gradient method. In [TH] the authors extended and 
improved the algorithm. The authors also applied the al- 
gorithm to the convex roof extension of the multipartite 
Meyer- Wallach measure [Tf] . 

In this paper we present an algorithm for the geometric 
measure of entanglement. The algorithm is easy to imple- 
ment, since every step is either the solution of an eigen- 
problem or finding a singular value decomposition of a 
matrix, and no further numerical techniques are needed. 

This paper is organized as follows. In Section |TT] we 
present the algorithm for pure and mixed states. We also 
discuss its properties and convergence. In Section [lll| we 
test our algorithm on bipartite and multipartite mixed 
states with known value of the geometric measure of en- 
tanglement. Further we compute an approximation of 
the geometric measure of entanglement for the isotropic 
states of three qubits, and the three-qubit XX model with 
a constant magnetic field. We conclude in Section |TV] 



II. ALGORITHM 



2. Multipartite states 

If we consider pure states \ip) on an n-partite Hilbert 
space H G ®^^iHi with n > 2, the geometric measure 
of entanglement is only known for a few special cases 
[IH] . In [ini [2n] the authors presented an algorithm for 
an approximation of Eq for pure states. For simplicity 
we will discuss the algorithm from [121110] for a pure state 
of three qubits, a generalization to arbitrary systems is 
done in the end of this section. 

Let It/i) be the given state of three qubits. The al- 
gorithm starts with a random product state |(/)o) = 
|0g^'')|0g^'')|0g'^'') of thrcc qubits, where the lower in- 
dex will be used for counting the steps of the algo- 
rithm and the upper index denotes the "number" of 

the qubit. Now we consider \^) = ((o[,^^|(o[)^^|) 

which is a pure unnormalized state on the space of 
the first qubit. If we want to maximize the overlap 
\{4>Q\%li)\ for fixed states jOg ) and |0q ), we have to 
replace loj,^^) by the state 10^^^) = , \ , U). The 
procedure is repeated for the second qubit, starting in 
the product state |0^^'') jOg^'') IOq'^') and resulting in the 

state |0^^^)|0f^)|0[,^^). Finally, the same maximization 
is done for the third qubit with the final state = 

|0i'')|0f')|0l^'>. In the same way we define the product 

state = |0i^^)|0l^^)|0i^^) to be the result of n iter- 
ations of the algorithm. In the following we will prove 
some properties of the algorithm. 

Proposition 1. Let |000) = lim„_j.oo he the product 
state after an infinite number of steps of the algorithm, 
then holds: 



Before we present our algorithm for general multipar- 
tite states, we begin with bipartite and multipartite pure 
states. 



A. Pure states 

1. Bipartite states 

For bipartite pure states G T-Li ® H2 the geometric 
measure of entanglement is given by pU| 



(lOOlV') = (OlOlV') = (001|V^) = 0. 



(12) 



EG{m)^l^Xl 



(11) 



where Xmax is the largest Schmidt coefficient of 1-0). Note 
that A^jjj. is also the maximal eigenvalue of Tri [|?/')(?/'|] 
and Tr2 Further let e Hi and \(j)2) & V.2 be 

the eigenstate corresponding to the maximal eigenvalue 
of Tr2 and Tri [!■!/)) (t/)!] respectively. Then the 

state 10) = |0i) ® 102) is a closest separable state to 



Proof. If (lOOj'f/') 7^ 0, then there exists a product state 
of the form |0) = |0(i))|OO) such that |(0|0)| > |(OOO|0)|. 
This means that |000) ^ lim„_>oo |0n), which is a con- 
tradiction to the definition of |000). Using the same ar- 
gument it can be seen that (0101-0) — (00110)) — also 
holds. □ 



From Proposition [T] we see that the state can be writ- 
ten as follows: 

It/-) = AilOOO) + A2IIIO) + A3IIOI) + A4IOII) -I- A5IIII), 

(13) 

where four of the coefficients A, can be chosen real and 



nonnegative, and X^i l-^i 



1. The form (^13 1 is also 



known as generalized Schmidt decomposition 

Proposition 2. The algorithm computes a generalized 
Schmidt decomposition of a pure state with an arbitrary 
given precision. 



3 



Proof. In order to find a generalized Schmidt decomposi- 
tion with a given precision e we need to find five param- 
eters with X]i=i IMjI^ = 1 a product basis {\ijk)} 
such that the state 



IV'approx) = Ml|000)+/i2|110)+Ai3|101)+Ai4|011)+Ai5|lll) 

(14) 

is closer to \^p) than e, i.e. -D dV") jV'approx)) < £ with 

the trace distance D {\^p) , {(f))) — y^l — KV'I^)]^. This is 
accomplished by the state 



(15) 



where |ijfc)„ = |*n^'')|jr?'')|fci^'') are the basis states 
after n iterations of the algorithm. The coeffi- 



Since in the present stage of the algorithm we only have 
the knowledge about the state |000) = |0(i))|0(2))|0(3)), 
the state can be computed as follows. Starting 
from the state \ip) we compute the unnormalized state 
I a) = \tp) — 1 000) (0001-0). In the second step we compute 
1/3) = |a) - E»<i |0W0(J'))(0W0(^'V)- In the final step 
we get iVi) = \P) - |0(^))(0(*)|/3). The state l^jj) is an 
unnormalized pure state of three qubits, and according 
to Proposition [l] applying the algorithm to it will give 
us the desired product basis {|ijfc)} with the property 
(211|'0) = {I2l\ij) = (11211/;) = 0. The expansion of the 
state in the final product basis {\ijk)} is a gener- 
alized Schmidt decomposition of \ip) [22]. Let {\ijk)n} 
be the computed product basis after n iterations of the 
algorithm. The approximated generalized Schmidt de- 
composition of \4>) becomes 



cients bijk are defined as follows: b 



100 



'010 



\ijk)n 



(18) 



with = biji = bjii = for i < j and 6^^ = a^-fc 
otherwise. N assures normalization of IV'n)- The preci- 
sion of the approximation is then given by D \ipn)) ~ 



In the same 



way we can find a generalized Schmidt decomposition for 
any multipartite pure state with an arbitrary precision. 



B. Mixed states 

The main idea of the algorithm for mixed states is a 
consequence of the fact, that the geometric measure of 
entanglement may also be written as j 23] 



Eg (p) = 1 - max F (p, a) 



(19) 



booi — 0, and bijk — ((V'Nj^)")* otherwise. 
N assures normalization of iV'n)- The trace dis- 
tance between lip) and Itf^n) becomes D{\^),\ipn)) = 

^|(V'|100)„|' -t- |(V|010)„|' + |(0|OO1)„|'. Using Propo- 
sition [1] we see that lim„^oo -D dV"), IV'n)) — 0. The 
wanted approximation IV'approx) is obtained by a state 
\iPn) such that D (IV-), <£■ □ 

Thus we showed, that the algorithm presented in the be- 
ginning of this section computes a generalized Schmidt 
decomposition of the given pure state. As the generalized 
Schmidt decomposition is in general not unique 2T1 122j , 
the result of the computation may depend on the choice 
of the initial product state \4>o). In particular, the final 
overlap 1 — [(OOOjV;)! does not have to be the geometric 
measure of entanglement, even for an infinite number of 
iterations. 

Finally we note that all results presented in this sec- 
tion can be extended to an arbitrary number of qubits. 
Then the equations have to be changed accordingly. For 
four qubits, eq. ^ becomes (lOOOj^A) = (OlOojiA) 
(OOIOI?/;) = (OOOll^f/^ 0. Moreover the results even hold 
if the subsystems are not qubits, but have arbitrary di- 
mensions. For simplicity we consider a pure state of three 
qutrits in the following. Again |000) — lim„_j.oo \4'n) de- 
notes the product state which is achieved after infinite 
number of iterations. Using the same arguments as in 
the proof of Proposition [T] we see: 

(lOOlV-) = (OIOIV-) = (001|?A) = 0, (16) 
(200|^) = (0201^;) = (0021?/^) = 0, (17) 

where |1) and |2) are arbitrary states orthogonal to |0) 
on the corresponding subspace. In order to find a gener- 
alized Schmidt decomposition we also have to find spe- 
cific states |1) and |2) for each subspace. Let \'ip) — 
X]i=o Sj=o X]j=o '^ijk\ijk) be the expansion of the state ^ 

in a product basis containing |000). Then consider the with pure separable states G 5, probabilities q^, a 



where S denotes the set of separable states and F {p, a) = 

(Tr [v^VpcrVp])^ is the fidehty. Let eH(g) Ha be a 
purification of p. It can be written as 



(20) 



with probabilities pi and p = J2iPi\'^i){'^i\- According 
to Uhlmann's theorem [3] page 410] and using (19) we 
can also write 

Ecip)^!-^^ max \{m\\ (21) 

where the maximization is done over all states \(f>) S 'H(^ 
Ha which are purifications of a separable state. Note 
that any \(f)) can be written in the form 



(22) 



unnormalized state 



J2i=i J2j=i J2j=i 0'ijk\ijk). unitary U acting on the Hilbert space Ha, and 
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From (21 1 we see, that we can get an approximation of 



Eg by maximizing the overlap |('0|(?!')| over all states \(j)) 
of the form ( 22 1 . Our approach for this maximization is 
the following: 



1. For fixed and we find a unitary U in (22) 
such that the overlap K-i/jjc/))! is maximal. 



2. For fixed U and we find states \(f>i) in (22 I such 



that the overlap KV'I'^)! is maximal. Note that this 
is in general only possible for bipartite states. For 
multipartite states we compute \(f>i) such that the 
overlap does not decrease. 



3. For fixed U and \<j)i) we find probabilities Qi in (22) 
such that the overlap |(?/)|(/))| is maximal. 

Steps 1-3 are iterated until the increase of the over- 
lap I (^Z" !<?!>) I is smaller than a small parameter e > 0. 
When the algorithm stops, the approximation of the ge- 
ometric measure of entanglement is given by Eq ( p) — 
(i/ilf/)) , where |(/)) is the final state of the form (22 1. 
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In the following section we will discuss the properties of 
the algorithm. Note that the order of the steps presented 
above can also be changed without changing these prop- 
erties. 



However, there is a way to circumvent this problem as fol- 
lows. We apply the algorithm described in Section [ll A 2 
to the state with the initial product state \(f)i), thus 
getting a final product state The state \cf)[) is not 

necessarily the closest separable state to |'0^), however it 
will be closer to jV'i) than the initial product state \4'i)- 
But then, if we replace by \(t>\), we get a better ap- 
proximation of the geometric measure of entanglement. 
This can be seen by noting that for the overlaps of the 
purifications holds: |('i/'|</'')l ^ where in |(/)') all 

product states were replaced by |(/>-). 

In the last step of the iteration we fix U which was found 
in the first step, and the separable states which were 
found in the second step. Using the method of Lagrange 
multipliers we find the optimal probabilities: 



(25) 



Let En (p) be the approximation of the geometric mea- 
sure of entanglement after n iterations of the algorithm. 
We will now prove the main property of the algorithm. 

Proposition 3. The approximated value of the geomet- 
ric measure of entanglement never increases in a step of 
the iteration: 



C. Properties 

In the following we will discuss some properties of the 
algorithm presented above. In the first step the prob- 
abilities qi and the separable pure states \(t>i) are fixed. 
The product Kt/iI^)! can be maximized using Uhlmann's 
theorem [21 page 410], it is maximal if U is chosen such 
that holds: 



En+l (p) < En (p) . 



(26) 



\/Z4t[/t 



(23) 



where A is a matrix defined as A 
J2i j \/Pi1j{'Pj\'4^i)\'i') Note that equation (23) is 

the polar decomposition of A, which can be computed 
efficiently for any matrix A [24^ . 

In the second step of the algorithm we fix [/, which was 
found in the step before. The probabilities qi are also 
unchanged. In order to maximize the overlap |('0|0)| the 
separable states have to be changed to the states 
for which holds: 



m^^) = ^JpAm), (24) 

with the states jV'i) — '^^j=vJ2j''^ij\/Pj\i^j)i where Uij — 

{i\U\j) are elements of U in the computational basis, and 
p'^ > is chosen such that is normalized. For bipar- 
tite states this step is evaluated according to the 

If \ip'^ is multipartite, the 
cannot be found in general. 



discussion in Section II A 1 
closest separable state 



Proof. It is sufficient to show that the overlap of the pu- 
rifications KV'I*/')! does not decrease in any step of the 
algorithm. This is seen directly from the definition of 
the algorithm in Section |IIB| □ 



D. Implementation 

First we set a small parameter e > 0. The algorithm 

starts with a random decomposition {pi, |'!/'i)}i=i into cP 

elements of the state p = '^i=iPi\''Pi) {'4'i\ and a separable 

decomposition {qi, \(f>i)}i^i of a random separable state 

a = X]i=i 'li\4'i){4'i\j where we demand that pi > and 
qi > for all 1 < i < d^. The steps 1-3 from the Section 
|IIB| can be implemented as follows: 

1 . Find the singular value decomposition of the matrix 

A - E^,JVm{^J\i^^)\^){J\, i-e- A - VDW with 
unitary matrices V, W and diagonal nonnegative 
matrix D. Define U = , noting that (23]) is 

fulfilled. 



2. Define unnormalized states 



\ai) = X^^iiv^l^J 

3 = 1 



(27) 



5 



with Uij = {i\U\j). Compute p'^ = {ai\ai) and 
1?/)^) = —y=\ai) for all i. For bipartite states 



compute separable pure states |(/>^) G S such that 

(V'il^'j) = \/Ps (IV'i)) ■ For multipartite states find 
product states \ip'j) which are closer to ji/'i) than the 
states {(pi) computed in the step before. This can 
be done applying the algorithm presented in Sec- 
tion 



II A 2 to the state jV',-) with the initial product 



state 



3. Compute q[ 



p'k\W'i''k' 



After performing steps 1-3 define a new separable state 
a' = J2i Qi\4>i){4'i\, which is an approximation of the clos- 
est separable state to p. If F{p,a') — F{p,a) > e, set 
= IV'^), If/-*) = Pi = K and q, = ql for all i 

and go back to step 1, otherwise stop. The computed 
approximation is Eq (p) = ^ — F {p, f ') ■ 



E. Convergence 

One of the most important questions regarding algo- 
rithms computing entanglement is whether or not the 
algorithm converges to the exact value of the entangle- 
ment measure, at least for infinite number of steps. For 
a general multipartite state with more than two parties 
the algorithm will converge to the wrong value with some 
nonzero probability, depending on the initial separable 
state. This is due to the fact, that the algorithm for pure 
multipartite states presented in Section |II A2| does not 
necessarily compute the correct value [20^ 2^ ■ 

For bipartite mixed states there is no full answer to this 
question, and testing the algorithm on bipartite states 
with known geometric measure of entanglement we did 
not observe convergence to a wrong value. However it can 
be shown that for some states and some special choice 
of the purifications and the algorithm does not 
compute the correct value even after an infinite number 
of iterations. To see this we consider a separable state 
p € S with rank r such that any separable decomposi- 
tion of p has more elements than r. The existence of 
such states is assured [5]. Let now {pi, \4'i)Yi^i be a 
decomposition of p which is optimal among all decom- 
positions with r elements, i.e. the average entanglement 
Ei=i Pi^G (IV'i)) is minimal among all decompositions 
into r elements. Further let be the closest separable 

state to IV'i) and we also choose qt = P'!^^'!*^'^! 

l^^Pk\{i'k\<Pk)\ 

Now we start the algorithm with the decompositions 
{Pi,\'^i)Vi=i and {qi,\4>i)}i^i, as described in the pre- 
vious section. Then the unitary U which maximizes the 
overlap of the purifications \ip) — ^ ■ y/pi\i^i} ^ N) and 
l<A) = Y.jV^\<Pj) <» U^j) is given hy U = 1. In the 
second step the algorithm will maximize the overlaps 
{(pilipi), which are already optimal. The same is true for 
the last step of the algorithm, where the probabilities qj 



are optimized. Thus the algorithm preserves the initial 
separable state, and does not compute the correct value 
even for infinite number of steps. 

To avoid the problem mentioned above the algorithm 
should always start with a separable state chosen at ran- 
dom, i.e. with random initial probabilities qi and random 
separable pure states \4>i). Moreover, the number of ini- 
tial nonzero probabilities qi should be at least (dimV.) . 

In the following section we will test the algorithm and 
present some applications for states with unknown geo- 
metric measure of entanglement. 



III. APPLICATIONS 
A. Testing the algorithm 

1. Two qubits 

If p is a two-qubit state, the geometric measure of entan- 
glement is given by (|9|. We applied our algorithm with 
£ — 10^^^ to lO'^ random states of two qubits and tested 
the computed value Eq against the exact value given in 
([9]). The maximal deviation Eq — Eq from the exact 
value was 6 • 10~^^. The average number of steps made 
by the algorithm was 291. 



2. Isotropic states 



We also tested our algorithm on the isotropic states in 
dimension d x d, these are states of the form 



p = p|$+)($^ 



1^1 



(28) 



with the maximally entangled state j^"*") = ^ EiLo 1^*)- 
For these states an exact expression for the geometric 
measure of entanglement was given in [S], the states are 
entangled if and only if p > j^-i. We applied our algo- 



rithm to the state ( 28 1 for 2 < d < 3 with the parameter 
e = IQ-^^ for p = O.Oln and < n < 99. The difference 
between the approximated value Eq and the exact value 
Eo was always less than 10~^°. 

In order to do the test for d = 4 within a reasonable 
time some modifications had to be applied. First, we 
minimized only over decompositions into d^ = 16 instead 
of d"* — 256 pure states. Further, for d = 4 the test was 
done on entangled states only, i.e. for p — O.Oln with 
20 < n < 99. The difference between the approximation 
Eg and the exact value Eg never exceeded 10"^'^. The 
results are summarized in Table |ll There N denotes the 
average number of steps made by the algorithm. 
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d 


2 


3 


4 


Eg - Eg 




< 10-"^ 


< 10"^^ 


N 


80 


516 


2259 



Table I. Precision of the approximation Eg — Eg an d the 



average number of steps A'' for the isotropic states ( 28 1 with 
parameter e = 10~^^. 



For the cases tested above the algorithm always con- 
verged into the correct value of Eq within the precision 
given m Table |I] with a single run of the algorithm. Note 
that in general more than one run with different initial 
parameters should be done to avoid convergence into a 
wrong value. Further we see from Table |T] that the pa- 
rameter e should not be used directly to quantify the 
precision of the approximation, although the deviation 
from the exact value is very small. 



3. Four qubits 

In [26] the authors computed the geometric measure of 
entanglement for a class of mixed states of four qubits. 
We tested our algorithm on the state p {t) , which for t = 
is defined as the four-qubit cluster state 

ICL4) = ^ (10000) + 10011) + |1100) - |1111)) . (29) 

For t > the diagonal terms of p are left invariant, and 
the off-diagonal components decay exponentially with t, 
i.e. 



Again for t > the off-diagonal components decay as 
in (30 1. The test was done with t = O.Oln for all 1 < 
n < 100, the difference Eq — Eq was always smaller than 
10~^^. The results are summarized in Table |lll There 
N denotes the average number of iterations made by the 
algorithm. 



P(0) 


ICL4) 


\Wa) 


1^4) 


Eg -Eg 


< 10-" 


< 10-" 


< 10-^^ 


N 


12 


173 


126 



Table II. Precision of the approximation Eg — Eg and the 
average number of steps iV for the four-qubit states presented 
in the text with parameter e = 10-^^. 



Note that the optimizations above were done over pure 
state decompositions into 2"* elements instead of 2*. This 
reduction was needed in order to do the computation 
within a reasonable time. Moreover we note that for 
very small parameter t = 0.01 we sometimes observed 
convergence into a wrong value. This is due to the fact 
that for small t the state p{t) is almost pure. As was 
mentioned in Section II E the algorithm can converge to 
wrong values for pure multipartite states. In these cases 
the algorithm was started again with random initial pa- 
rameters. To get an impression we mention that for the 
last example p(O.Ol) the algorithm sometimes converged 
to Eg- Eg kS- IQ'^. 



4- Comparison with other algorithms 



Pkl (t) 



Pkl (0) 
e^'pkl (0) 



for k — I, 
for k^l. 



(30) 



-15 



We applied our algorithm with parameter e = 10 on 
the states p{t) with t = O.Oln for all 1 < n < 100. 
The discrepancy between the approximated value and the 
exact value given in [2S] was always smaller than lO-^''. 

The same test was done for the state p (t) , which for t = 
is defined as the four-qubit W state 

\Wi) = ^ (10001) + 10010) + 10100) -t- |1000)) , 

and for i > the off-diagonal components decay expo- 
nentially as given in (30 1. There the discrepancy be- 



tween the approximation and the exact value was always 
smaller than 10"^^. 

Finally we tested our algorithm on the four-qubit state 
p (t) , which for i = is defined as the symmetrized Dicke 
state 

\D4) = ^ (10011) + 10101) -I- |1001) -I- 11100) -I- 10110) 
v6 



11010)) 



(31) 



A significant difference between our algorithm and the 
algorithms presented in [121 US] is the fact, that our algo- 
rithm implies only the solution of the eigenproblem and 
finding a singular value decomposition. For both prob- 
lems efficient numerical algorithms exist [23], implying 
that each step of our algorithm can be done efficiently. 
The algorithms based on conjugate gradients usually im- 
ply a line search |15| . It is not known to us whether a line 
search can in general be done efficiently for the problem 
considered here. 



As noted in Section |III A 1| the average number of itera- 
tions made by our algorithm for random two-qubit states 
with parameter e — 10^^^ was 291. This is comparable 
to the performance of the conjugate gradient algorithm, 
for comparison see Figure 1 in [1^. 



B. On additivity of entanglement 

A measure of entanglement E is called additive, if for any 
two states p"^^ and a^^ holds [B]: 



E{p 



AB ^ ^AB-^ 



E(p 



AB\ 



E a 



AB\ 



(32) 
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where the entanglement between the parties A and B is 
considered. 

For pure states IV''^"^) and 10"^'^) we see that 

(1^^^) ® = (1^^^)) i^, (10-^^)) , (33) 

with Fs (p) ~ mSLxF (p, cr) and the fidehty F {p, a) — 

(Tr [\/y/P'^y/p] ) ■ From (33) we see that the geomet- 



ric measure of entanglement is not additive. Note, that 
for the entanglement of formation nonadditivity has also 
been proved [?7] . 

We consider the logarithmic entanglement 



Flog (p) = - log2 F, (p) 



(34) 



which is additive for pure bipartite states, as is seen 



from (331. In general holds: F, (p-^^ (g) cr^^) > 



Fs [p^^) Fs [a^^), and thus the logarithmic entangle- 
ment is subadditive: 

) < Flog [p 



T? I AB „ AB\ 

Flog [P ®(J 



0+i?zo,(0. (35) 



We use our algorithm to test the inequality ( 35 1 . 
Note that for two-qubit states p we get Fs (p) = 

We take p'^^ and a^^ to be 



random states of two qubits, and apply the algorithm 
to p"*-^ ® a^^ with parameter e = 10~^. This 
procedure is repeated 100 times, each time the com- 
puted approximation Fs (p^^ <8) a^^) was slightly below 



w.. . s yr^"" -^^^ 

3-4s~j ^^AB\^^ which means that we could not dis- 
prove additivity of logarithmic entanglement in this way. 



Fs 



The difference Fs (p^^) Fs 
always smaller than 10~^. 



(7 



AB\ 



Fs 



..AB 



) a 



AB 



was 




Figure 1 . Approximation of the geometric measure of entan- 
glement Ec for isotropic states of three-qubits given in ( 36 1 



as a function of p (solid line) compared to the two-qubit case 
(dashed line). 



with I GHZ) 



1 

V2 



(|000) + |111)). They are known to 



be fully separable if and only if p < | [IH]- We apply 
our algorithm to these states with parameter e = 10~^ 
for p > |. The result is shown in Figure [l] (solid line). 
The plot can be compared to the geometric measure of 
entanglement of the isotropic states of two qubits, see 
dashed line in Figure [l] In the limit p 1 the state 
becomes the pure GHZ state with Eq (|GHZ)) = 10. 



2. XX model 



C. Applications to 3 qubits 

In this section we apply our algorithm to 3-qubit states 
with unknown value of Fq. If d is the dimension of the 
total Hilbert space, then for any p there always exists an 
optimal decomposition with at most cP elements [25] . A 
decomposition {pi,\4'i)} of a state p = J2iPi\'^i){^i\ 
called optimal if its average entanglement is equal to the 
geometric measure of entanglement: X^i-Pi^G (IV'i)) = 
Eg (p). In order to make sure that the algorithm al- 
ways has the chance to find the optimal decomposition, 
all minimizations in this section were done over decom- 
positions into (P = 2^ = 64 pure states. In order to do 
the computation within a reasonable time we used the 
parameter e ~ 10^^. 



1. Isotropic states 
Isotropic states of three qubits have the form 



p|GHZ)(GHZ| 



-1, 



(36) 



As a final example we apply our algorithm to the isotropic 
XX model of 3 qubits in a constant magnetic field. The 
corresponding Hamiltonian is given by |29L I30| 



H 



B 



2 ^ 

i=l 



J 



3 

E 

1=1 



cr, (T, 



i+l 



i+lj 



(37) 



with periodic boundary conditions: erf = af and 
af. In thermal equilibrium the system is found in the 



mixed state p — with Z = Tr e"*^ . In the fol- 

lowing we set k ~ I. 

The results of the approximation with parameter e = 
10~^ are shown in Figure |2] They can be compared to 
the results for two qubits in Figure 4]. For different 
values of the magnetic field B we observe a different be- 
havior of the system in the low temperature limit. This 
behavior will be explained in the following. 

Note that the Hamiltonian ([37| has four nondegenerate 



eigenvalues ±f-B, and 4 J ±^5. Further the following 
two eigenvalues are degenerated twice: — 2 J ± ^B. For 
vanishing magnetic field the ground state of the system 
is a mixture of the four eigenstates corresponding to the 
eigenvalue —2J with equal probabilities. In this case we 
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Figure 2. Approximation of the geometric measure of en- 
tanglement Eg plotted as fu nct ion of the temperature T for 

The parameter J is set to 



P = 



kT 



and k 



with H given in (37 1 
1. 



responding to the eigenvalue — 2J — ^B. As can be seen 
from the dashed curve in Figure [2] for T — > the approxi- 
mation becomes Eq ~ | ™ this case. In the case B — 2J, 
there are three eigenstates corresponding to the smallest 
eigenvalue — 3J. The approximated value for T — > in 
this case becomes Eg « 0.116, see dotted curve in Figure 
|2] Finally, for B > 2J the ground state is the product 
state I 111), and the entanglement vanishes for T ^ 0, as 
is seen from the dot-dashed curve in Figure [2] 

In Figure [s] we show the plot of Eq as function of the 
magnetic field B for three different temperatures T. For 
T — > we observe that Eq becomes a nonanalytic func- 
tion of B for two different values of the magnetic field, 
namely for i? = and B — 2J. This is a significant differ- 
ence to the two-qubit case, where such behavior occurred 
only for a single value of B [311, Figure 5] . 



IV. CONCLUDING REMARKS 



Eg 
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Figure 3. Approximation of the geometric measure of entan- 
glement Eg for fixed values of T plotted as function of the 
magnetic field B. The parameter J is set to |, and k = 1. 



get Eg ~ J for T — 0, see solid curve in Figure [2j For 
small nonzero magnetic field < B < 2J the ground 
state of the system is the mixture of the eigenstates cor- 



In this paper we presented an algorithm for approximat- 
ing the geometric measure of entanglement for arbitrary 
multipartite mixed states. The algorithm is based on a 
connection between the geometric measure of entangle- 
ment and the fidelity fl^. It is easy implementable, since 
it implies only the solution of an eigenproblem and find- 
ing a singular value decomposition. We tested our algo- 
rithm on bipartite and multipartite mixed states, where 
an exact formula for the geometric measure of entangle- 
ment is known. In all cases we found convergence to 
the exact value. For two qubits, the performance of our 
algorithm is comparable to the performance of the algo- 
rithms based on conjugate gradients. We also applied 
our algorithms to the isotropic state of three qubits, and 
the threc-qubit XX-model with external magnetic field. 

In our tests on bipartite mixed states with known value 
of the geometric measure of entanglement our algorithm 
always converged to the correct value within a given pre- 
cision. It remains an open question whether this is al- 
ways the case. For quantum states with more than two 
parties the algorithm can converge to wrong values with 
nonzero probability. In general more than one run of 
the algorithm with different initial parameters should be 
performed. 
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